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1 Introduction 



It is a typical phenomenon in thermal field theory that the conventional per- 
turbative loop expansion breaks down because of infrared singularities which 
are usually cured by the generation of thermal masses[l]. In the particularly 
simple case of a massless A0^-theory, the leading temperature contribution to 
the self-energy consists of a simple mass term with m = y^T. An improved 
perturbation theory requires that this thermal mass be treated on a par with 
the tree-level inverse propagator whenever the momentum of the latter be- 
comes comparable to (or smaller than) the former. A consistent way of doing 
so is to add the thermal mass term to the bare Lagrangian and to subtract it 
again through counterterms at higher loop orders. 

In general, however, the thermal self-energy is more complicated than a con- 
stant mass squared and it may also be the case that higher vertex functions 
are equally important. This is the case in QED and QCD, for which a re- 
summed perturbation theory has been developed by Braaten and Pisarski[2]. 
The leading temperature contributions to be resummed have been obtained 
by Frenkel and Taylor [3] and independently by the former, who have coined 
the term "hard thermal loops" . These hard thermal loops (HTL) turned out 
to satisfy ghost-free Ward identities[4,5]. Manifestly gauge invariant as well 
as rather compact expressions for the effective Lagrangians which summarize 
them were found in Refs. [6-9]. 

This resummation program has been applied successfully to a mimbcr of prob- 
lems (see e.g. the forthcoming review in Ref. [10]). Occasionally, however, it 
turned out that the thermal masses contained in the HTL effective action were 
insufficient to screen all infrared singularities. For example, through quasi- 
particle mass-shell singularities, the damping rates of plasma excitations with 
nonzero momentum as well as the next-to-leading order screening corrections 
are blown up by unscreened magnetostatic modes [11]. However, assuming 
the existence of some effective infrared cutoff at the (unfortunately entirely 
nonperturbative) magnetic mass scale at least allows one to obtain a leading 
logarithmic correction. 

Another shortcoming of the conventional HTL resummation has been revealed 
recently in the attempt to calculate the production rate of soft real photons 
from a quark-gluon plasma. In the case of hard real photons, a complete 
leading order calculation was carried out in Ref. [12], where it was found 
that the mass singularities of the quarks are screened by HTL corrections. 
But with soft real photons, the hard thermal loops themselves introduced 
uncancelled coUinear singularities [13]. We shall concentrate on this latter 
type of singularities in what follows. 
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In a quite different context a similar difficulty was encountered in Rcf. [14]. 
In the simpler gauge theory of scalar electrodynamics, the complete ncxt- 
to-leading order dispersion laws of photonic excitations could be obtained in 
a resummed one-loop calculation. However, the perturbative result for the 
longitudinal branch turned out to become unreliable as the light cone was ap- 
proached with increasing plasmon momentum. The reason for this breakdown 
was again a collinear singularity of a HTL diagram. There it was found that 
it can be removed by a further resummation. 

In this paper, we extend this strategy to the case of QED and QCD and 
we shall show that a manifestly gauge invariant HTL effective action can be 
found which is completely free from collinear singularities. It improves upon 
the effective action of Refs. [6-8] with which it coincides for external momenta 
that are sufficiently far from the light-cone. 

In the next section we argue that the hard propagators have to be dressed 
by asymptotic thermal masses whenever collinear singularities make the hard 
thermal loops themselves sensitive to such higher-order corrections. In sect. 3 
we carry out a resummation of the asymptotic thermal masses first for the 
case of purely gluonic QCD, and in sect. 4 when fermions are included. The 
analytic structures obtained differ, but in each case the improved hard thermal 
loops can be summarized by a gauge invariant effective action. Sect. 5 contains 
our conclusions and gives an outlook to potential applications. 



2 Hard thermal loops in the vicinity of the light-cone 

Consider the 00-component of the polarization tensor (i.e. self-energy) at lead- 
ing order as given by its hard thermal loop. It reads universally [15] 



with m = eT/3 and n is the Bose function, where e is either the coupling 
constant of QED or 




Qo - pq/p 



Qo 






(2) 



for SU(iV) with Nf fermions. 
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Fig. 1. Prototypical Schwinger- Dyson equation — for simplicity for a 
(^^-theory. Internal propagators are always dressed, whereas vertices ap- 
pear in bare and dressed versions. 



Close to the light-cone IIoo is found to diverge logarithmically, 

Hoo In - with = Q'^/q^ . (3) 

£ 

The origin of this divergence is a coUinear singularity in the loop integral (1) 
at pq = pq when Qq — q. 

The next-to-leading order contribution for soft momenta Q = (Qo, q), Qo, q ~ 
eT is given by resummed one-loop diagrams. In both scalar electrodynamics 
[14] and QCD [16], it turns out to diverge hke 

^IIoo ^rn?- (4) 



when approaching the light-cone. Evidently, the resummed perturbative series 
breaks down for £ < e, for then ^IIoo > Hoo- 

In the case of scalar electrodynamics (massless and without scalar self-interact- 
ions) it is relatively easy to analyse this problem [14] . The coUinear singularity 
in (3) is brought about by the masslessness of the internal scalar particles. The 
even stronger singularity in (4) is generated by a premature restriction to soft 
loop momenta. The scalar propagators are dressed and therefore massive, but 
one has to subtract off the contributions already covered by the bare HTL 
diagram. This again involves massless scalar propagators, which become large 
at the light-cone irrespective of the momentum scale. Thus (4) reveals a latent 
UV- divergence. 

However, from the Schwinger-Dyson equations for the various Green functions 
(Fig. 1) it is clear that in the full theory the scalar lines are always massive 
with mass squared = |e^T^(l + 0(e)). The coUinear singularities of Hoo 
at the light-cone, which are produced by scalar particles, should therefore be 
spurious. 

Indeed, keeping the thermal masses also for the hard scalar propagators yields 
a finite result for the HTL contribution to Hqo where e in (3) gets replaced by 
e. The increasingly singular contributions (3,4) are thus seen to get cancelled 
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by higher loop diagrams corresponding to a hard scalar loop where scalar 
self-energy diagrams are inserted repeatedly. Away from the light-cone, such 
insertions are suppressed by powers of e^, but they cease to be so when e < e. 

Resumming thermal masses already at the stage where the hard thermal loops 
are being calculated raises the question about the systematics of this proce- 
dure. For this it is useful to think in terms of a renormalization group ap- 
proach [17]. The conventional resummation program of Braaten and Pisarski 
can be understood as a two-step procedure. First an arbitrary scale A with 
T » A » eT is introduced which divides all momenta and energies in hard 
and soft, and an effective theory is built from integrating out all the hard 
modes. To leading order, the result is independent of the actual value of A, 
and is given by the HTL effective action. In a second step, Green functions 
with external soft momenta are calculated using this effective action with soft 
loop momenta. The degree to which a perturbative expansion based on the 
effective theory can possibly make sense is limited by the accurateness of the 
effective action. Higher loop orders of the former require sufficiently high loop 
orders of the latter. 

In the above example it became apparent that for soft momenta with s < e 
the effective action receives contributions from hard diagrams with arbitrarily 
high loop orders. In order to sum them systematically one would have to solve 
the Schwinger-Dyson equations, which are an infinite set of coupled equations. 
However, the present singularities are caused primarily by the masslessness of 
the bare propagators, whereas the structure of the Schwinger-Dyson equations 
is such that only full propagators enter (in contrast to vertex functions which 
appear in both bare and dressed forms). As long as it is sufficient to sum only 
self-energy insertions to produce massive propagators, there is no particular 
problem with overcounting. 

In fact, also in the more complicated gauge theories like spinor electrodynamics 
and QCD the physical degrees of freedom retain thermal masses for momenta 
p ^ eT. In the case of gauge bosons, there is a transverse and a longitudinal 
branch of quasi-particle mass-shells, but the longitudinal branch rapidly dies 
out with increasing momentum. In the HTL approximation, the longitudinal 
plasmons approach the light-cone, but in doing so the corresponding residue 
vanishes exponentially [18]. The transverse mode, on the other hand, tends 
toward an effective asymptotic mass 

ml ^ H,(Q^ = 0) = ^ + O(e^rA) (5) 

and residue 1. (Recall that in QCD e is defined by (2).) Likewise, the disper- 
sion laws for ultrarelativistic fermions have a collective branch (the "plasmino" 
with a flipped relation between helicity and chirality) which dies out with in- 
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creasing momentum and a normal one that remains, again with an asymptotic 
mass proportional to e^T^ + O(e^TA) (without spoiling chirality) [19]. 

The HTL effective action is generated entirely by the physical degrees of free- 
dom of the bare theory, which is particularly evident when the Coulomb gauge 
is used, where only transverse modes are heated. ^ To leading order, these con- 
tribute only as real particles living on their mass-shell, which is the light-cone 
= 0. Correspondingly, the coUinear singularities present in conventional 
HTL diagrams with external light-like momenta are all removed when the 
above asymptotic thermal masses for transverse gauge bosons and any other 
ultrarelativistic particles are included. 

It is an important feature of the asymptotic thermal masses that the result 
(5) as well as the corresponding one for fermions (see below) does not depend 
on the condition Qo,q <^T which is otherwise essential for the derivation of 
the HTL results. It holds for arbitrary momentum as shown in Eq. (A. 21) of 
Ref. [21]. Moreover, the definition (5) is itself not sensitive to a soft modifica- 
tion of the hard propagators, because H^ is not singular at the light-cone. 

However, the smallness of the asymptotic masses might help formerly negligi- 
ble contributions of higher-order diagrams to increase when approaching the 
light-cone. This may raise the necessity to reorganize the perturbation series. 

In gauge theories one might expect that giving masses to the transverse gauge 
boson modes without modifying also their vertices might violate gauge invari- 
ance. Nevertheless, we shall see presently that gauge invariance of the HTL 
effective action is maintained if the above procedure is followed to generalize 
hard thermal loops to soft hghthke momenta. 



3 Purely gluonic QCD 



In this section we shall first treat purely gluonic QCD, where = g'^N. 
The improved HTL diagrams are obtained almost exactly as usually, the only 
difference being that the gluon propagator for hard momenta p > A is modified 
to take into account the asymptotic thermal mass for transverse gluons, 

GAP) = 

+ {B^, + aL*^.) Ao (6) 



• 1 A 1 A 1 



oo 



^ By employing a nonstandard real-time formalism, this feature can be kept also 
for covariant gauges [20]. 



6 




Fig. 2. The contributions to the improved gluon self-energy. The wavy 
hnes with a blob represent gluon propagators whose transverse piece 
resums the asymptotic thermal mass; the dotted line represents the (un- 
changed) ghost propagator. 



and — — 3^,^ — Z)^,^ , 5^,^ = , D^iv — ~p2~ 

where = Uf^P"^ — Pn{PU) and [/^ = (1, 0) is the four velocity of the thermal 
bath at rest. 



3.1 Gluon self-energy 



Since the collincar singularities in the bare one-loop polarization tensor are 
only logarithmic, the extraction of the contribution with soft external mo- 
menta can proceed as usually (see, however, the Appendix for a pitfall). Soft 
external momenta in the numerator can be neglected when compared to hard 
loop momenta, and the familiar expression is obtained, but with in place 
of the originally massless propagators. 

Let us demonstrate this in the example of IIqq{Q). With the modified gluon 
propagator (6), the diagrams shown in Fig. 2 respectively yield 

nor(Q) - g'N^^ {^m - 4A- A^Poi^o - Ao AoPoi^o} 
IllfiQ)^-g'N^^3Am (8) 
nt^\Q) = g'N^^A^AoPoKo, 

such that in the sum only massive propagators appear, 

noo(<5) = g'N^^ {-2A^ - 4A- A^PoXo} ■ (9) 



In the following we use the imaginary-time formalism, so the zero components 



7 



of the momenta P = {Pq, p) are discrete Matsubara frequencies Pq — 2mnT. 
The symbol f,p is defined as 



Throughout our paper Q is the external momentum, P is summed over and K 
is the difference K — Q — P. An index ~ means the transformation P ^ K, 
e.g. Ao = 1/K\ 

Performing the sum over the Matsubara frequencies yields —2g^Nfp A^^ — 
with = and 



Tp Tp \ Pn-Ur, Po + UJ„\ I An - U)k 



'P [Po-Up Po + UJpj IKq- Uk Kq + UJk . 

d^p / 1 1 \ 



/dp 
[n{ujp) + n{ujk)] 



(27r)3 XQo + ujp + ujk Qo- cOp- UJk J 
+ [nM - nM] ( ^ - ^- ^— - ) (11) 

with 0;^ = + m^. After the frequency sum is done, we continue Qq to real 
values, keeping a small imaginary part +is in mind. 

Since cUp — a;^ is a soft quantity of order q, 

u;p-u;k-zq-zq^-f{l-z^) + o('-^], z^^, (12) 
2p^ 2p \P J PQ 



it is useful to expand n{LJk) = n{up) + n'{up){uk — uip) + 0{q^). Moreover, 
the first two denominators on the r.h.s. of (11) are far from potential zeros 
for hard loop momenta, where Qq is negligible against uip + uik, so one can 
approximate cjp, Uk ^ p there. The remaining denominators are those that 
give rise to collinear singularities at the light-cone when moo is neglected. 
However, since these singularities are only logarithmic, we may also expand 
their numerators, dropping the terms that have too little power at large p to 
contribute to the T^-part. This leads to 

Hoo = Sm' + 2g'N [ ^ n'{p) ——^ . (13) 

J (27r)^ - [Up-UkY 

Sufficiently far from Q'^ = 0, one can replace ujp — ui^ zq. In the vicinity 
of the light-cone, however, the main contribution comes from \z\ ~ 1 and 
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the second term in (12) becomes important, whereas the subsequent term is 
suppressed by (1 — 2;^). This latter term is effectively of 0{q^/p^), because 
the dominant contribution is generated when (1 — z^) ~ w? j]? . We therefore 
obtain 



noo(Q) = 3m^ + / ^ n'(p) ^ ^ (14) 

— 3m^ — %; (da-. -r- - f dz , ^ - (15) 

3 (e"-l)^2y 0^g„ _ 

where /x^ = !^ = 

Although the angular integration is no longer decoupled from the integral over 
the modulus of the loop momentum, this is still very close to the structure 
found in the conventional hard thermal loop, see Eq. (1). The difference is 
only that Q/^ does not appear in the scalar product with a light-like "unit 
vector" which is integrated over, but is contracted with a timehke one whose 
zero component is subject to a certain averaging. 



Introducing 



Y^{a) = {Yo{a),e) with Yo{a) = (16) 



a 



and the double averaging 



3 7 a^gc I 







allows us to write the improved HTL polarization tensor in the same compact 
way as the conventional one. 



n^, = 3m2 (U^U, - ^Y^Y, \ . (18) 



In this expression Yq ^ 1 is essential only in the denominator, since n g. 
Strictly speaking, the lower bound on the integration variable a is given by 
A/T, but this is negligible when concentrating on the leading contributions 

is now defined also for soft hghthke momenta. The former singularity of 
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rioo at = has disappeared and is replaced by ln(const./g'), to wit 



g^NT^ 2 , 1 _ , C(2) 1 



noo(g^ = 0) = -^^<|ln- + --7+^| . (19) 



There is still a logarithmic branch cut for \Qq\ < q. The imaginary part along 
this cut, however, is a smooth albeit nonanalytic function in the momentum 
variables, 

^mUooiQ)^^^ 1 da . f^" ^ forg2<0 (20) 

ao{Q) 



and with ao{Q) = fJ^Qo/ V~Q^' ■ \Q^\/Ql 3> g, the lower integration bound 
ao{Q) is of order g, hence neghgible at leading order. However, it becomes of 
order 1 when \Q'^\/Qo ^ g, and finally goes up to infinity for \Q'^\ — > 0, thereby 
bringing (20) down to zero. 



3.2 Gluon vertex functions 



For gluon vertex functions more complicated analytic structures arise. Per- 
forming the analogous steps as above, we find e.g. for the 000-component of 
the 3-gluon- vertex function 

rooo(Q) -R) —Q — -R) 

1 1 



+ 



Qo + '^p-q ~ '^p -Rq + ^p-q-r ~ ^p-q 
1 \ 



Qo ~ <^p-q + ^p Ro ~ ^p-q-r + ^p-q, 
1 



~{Qo + Rq 



^p—q—r ^p 



Qo — LOp^q + LOpQo + Ro — LOp^q^r + ^pjj ^ 

after dropping all contributions <^ g^T^- 

Without the asymptotic thermal masses there are coUinear singularities when 
any of the external momenta is light-like. These are again only logarithmic 
as long as the external momenta have different directions. Assuming this, one 
can for each denominator in turn expand the energy differences according to 

2 

(12). This again amounts to modifying Qq — > Qq(1 + ^) and likewise Rq. 
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After that the two pairs of terms in the round brackets can be combined by 
changing p — > — p. 

As a result, the modified gluon vertex can be cast into the same form as was 
possible for the conventional HTL vertex 

r™ = 3gm I - (7(0^); . (22) 

but with the averaging and the vector Y redefined according to (17) and (16). 
Evaluating the angular integral gives 

Tooo = Sgm^ I \RoM {Q,R)-{Qo + Ro)M{Q,Q + R)] (23) 

where 



Q = (go^/l + /^V«^q) (24) 



and M. the Lorentz- invariant function introduced in Ref. [3], 



In the conventional HTL result there are logarithmic singularities at = 0, 
= 0, and (Q + i?)^ as well as (generally nonsingular) branch points at 
A{Q,R) = A{Q,Q + R) = 0. In the improved result (22), these singularities 
are removed because in the a-average they are smeared according to (24). 

As concerns higher vertex functions, the above argument that led to the form 
(22) can be essentially repeated. 



3.3 Improved effective action 



Because of the formal similarity of the improved hard thermal loops to the 
conventional ones — we only had to redefine the 4- vector that is used in the 
angular average and extend the averaging — it is natural to guess that one 
can also take over the compact effective action by Braaten and Pisarski [8] 

5eff = -3m'\ J d'xFrix) (^^) F,Mx) , (26) 
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but now with the definitions (16) and (17). These modifications obviously do 
not interfere with the manifest gauge invariance of the HTL effective action. 

One can readily verify that (26) contains the improved version of the HTL 
two-point function, but some manipulations are needed to bring it into the 
form derived above (see the Appendix). 

In the case of vertex functions, it is advantageous to use the earlier version 
of the HTL effective action of Taylor and Wong [6], which is, however, not 
manifestly gauge invariant. We therefore demonstrate that the redefinitions 
(16,17) do not interfere with its actual gauge invariance either. 

We first write 

oo 

S,s = T.^n[A], (27) 

n=2 

where Vn collects the contributions n-linear in A^. Prom the improved gluon 
self-energy we have 

S2[A] = 1 1 ^^"'^(-Q) n,.(g) A-(g) 

^Sm'' Jd^xtr (^[Ho{x)]^ - H{x)^H{x)j (28) 

with the abbreviations H{x) = Y^'A^^{x)T'' and Ho{x) = YoA^{x)T''. With 
the ansatz 

S = 3m^ Jd^xtr {{[Ho{x)]^) - • (29) 

one can find a gauge-invariant functional with "boundary condition" (28) 
as follows. 

Under infinitesimal gauge transformations with parameter lj we have 

Si^Hq — YoOqu; — ig[Ho, lo] and S^H — Ydu — ig[H, u] —: Vu . (30) 

Gauge invariance of S implies = for all uj. This requires 

2{Y^doH) = {VdH(t>[H]) . (31) 
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Up to terms that average to zero, a solution to this equation is given by a 
functional satisfying 

S„m\-iYo9oH wth i = ^_J-^J^. (32) 



Note that here there is a small difference to the conventional case. Because of 
the new averaging procedure the inverse of (Yd) always exists. 

Counting only explicit coupling constants (not those implicit in Yq which de- 
pends on II (X g) , we have gdg(f) = J d'^xtr HdH(f>- That is 

f °° r 1 1 " 1 

gdg<P[H] = 2jd'xiiY. ^XYQ^^'aH, *]| y^YodoH . (33) 



Integrating with respect to g yields 



S,s[A] = 2>m^ J d^xtr (^{Hof + {YodoH)F (^[^^, *]) (34) 



oo 



where F(z) = 2 ^ 

n=0 ^ + 2 



This functional is gauge invariant by construction. In fact, it is identical in 
form with that of Ref. [6], only the meaning of the symbols has changed. It 
reproduces the improved vertex functions in the form in which we had them 
obtained in the previous section, since Yo's in the numerators can be put to 
1. This just drops terms that are suppressed by powers of g, which we have 
always discarded. 

We thus have shown that inclusion of the asymptotic gluon mass (5) removes 
the coUinear singularities of the hard thermal loops without spoiling gauge 
invariance. 



4 Inclusion of fermions 



Defining the two structure functions of the fermion self-energy E at finite 
temperature by 

E = aQoTo + &q7 , a = tr 70E b = --^ tr q7E , (35) 

4:Qq 4^^ 
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the dressed propagator reads 



(1 + a)Qo7o - (1 + b)qj 

(l + afQl-il + byq^ ■ ^ ' 



The leading high-temperature contributions are 

„__^JQ^] _ ^^MlU^9iJQi±i]\ (37) 



2g(5o \Qo-qJ ' q"^ [ 2q \Qo-q, 



with Mj = ^g^CfT"^ (the QED case is covered by replacing g'^Cf e^). 
For soft momenta, this contains an extra collective mode, the "plasmino". 
However, as with the plasmon, its residue vanishes exponentially fast with 
increasing q. For Qo,q^ gT, the dressed propagator approaches 

Q2-Ml + 0iMyq^)' ^^^^ 



where 

= 2{Qla + q%) = 2Mj. (39) 



Like the asymptotic gluon mass, this latter result (in contrast to (37)) is 
even an exact one-loop result when = 0, i.e. not merely the leading high- 
temperature term for Qq, q <^T (see the appendix of Ref. [22]). 

With only external gluons the fermion loop produces the same HTL contri- 
butions as the purely gluonic one-loop diagrams did, except that instead of 
there is an overall factor of Nf/2. For lightlike or nearly lightlike external 
momenta, the asymptotic fermion mass M^o becomes important and appears 
just in place of moo- 

However, the situation is essentially different when hard thermal loops with 
external fermions are considered. The only hard thermal loops are those with 
two external fermion lines. Those are given by a loop which involves necessar- 
ily both internal fermion and gauge boson propagators, which have different 
asymptotic masses. 



4-1 Fermion self-energy 



Dressing the hard propagators in the fermion self-energy with their respective 
asymptotic masses gives the improved hard thermal loop 
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m) = --/c, j ^ Hp) + n{p)] (70 - ^ 

1 



X <^ + \ (40) 

— UJp + i^p-q Qo+ COp — UJp-q J 



where h is the Fermi- Dirac distribution function and ujk = yk"^ + 
In the difference 



u;,-Co, = zq+'^^^-^-f{l-z') + 0{^-^), z=^, (41) 
2p 2p p^ pq 



the dominant term that shifts the poles in the integrand of (40) is now inde- 
pendent of Q and proportional to the difference of the two asymptotic masses. 
The other correction in (41) that is of comparable magnitude is suppressed 
for \z\ ~ 1, which makes it effectively of order q^/p"^. 

The two terms in the curly brackets of (40) therefore do not combine when 
exchanging p — > — p and we arrive at an improved HTL fermion self-energy 
of the form 



E(Q) - -Mj I { + (yg)!:^ 



(42) 



where 



i; = F,(oo) = (l,e), l--.jd<^^^^ dm^ ^ - .(43) 

°' 



So the 4-vector Y is light-like as with the conventional hard thermal loops. 
The coUinear singularity is instead smeared out by the addition of idm/o; 
with a slightly modified averaging prescription for a. 

The singularity of KeS at the light-cone is again cut off by the asymptotic 
masses with the result 



»eaW-O) = -0{„>^ + |.n2 + l-, + £|l} (44) 

^eb(Q^ = Q) = ^ea(Q^ = Q) + ^ (45) 

Q 

However, in contrast to the purely gluonic case, the imaginary part of the 
fermion self-energy is left unchanged for all spacelike momenta and remains 
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nonvanishing at the light-cone, 



Qma{Q^ = 0) ^^mb{Q^ ^0) ^ (46) 

z q 



dropping rapidly but smoothly to zero only for timelike > g'^{gTy 



4-2 Improved effective action with fermions 



In the case of the conventional hard thermal loops, the effective action for 
fermionic Green functions is given by the simplest possibility to bring its 
bilinear part, which is determined by the fermion self-energy, into a gauge- 
invariant form. With equal ease, this can be done with the improved HTL 
fermion self-energy, 

£.» = -^i^M] I I + ,,^(.) (47) 

^ a a ' 



with D^'^dn- igT'^A'^i'. 

This obviously requires that not only the fermion self-energy but all higher 
vertex functions can be broken up in two parts which differ only by the sign 
in front of dm. This is indeed the case as can be shown by induction. 

Let us first explain how the structure of (42) arises. This simplest of the 
fermionic Green functions contains one fermionic and one bosonic propagator 
inside the loop. Decomposing each in partial fractions gives rise to a sum of 
products of the form 

Aj.Ab ^ y -—^ -—^ (48) 

^2 -H) + (yiOJf Bo + a2UJb 

where wc have denoted the respective momenta by and and the ai are 
signs. F and B differ only by soft momenta, so we choose their orientation such 
that Fn ^ and keep this convention when we shall be considering ad- 

ditional propagators in the loop. Performing now a further partial fractioning 
on (48) yields 



(T1CT2 ^' 



'<o- Fo + aiLVb - CF2UJf \Fq + aiLJf Bq + a2UJb^ 



Summing over the Matsubara frequencies will give —ain{u}f ) in the first, and 
—ain{u}b) in the second term. Of the four parts of the sum only two can 
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contribute to hard thermal loops, namely those where n{oOf) and n{u}b) acquire 
the same sign, di = (T2- Otherwise the contribution will be suppressed by a 
hard denominator in the prefactor. 

Now let us consider a loop with one additional propagator, Api or A^/. This 
brings in another factor 1/{Fq + cr'uf) or 1/{Bq + a'uy). Decomposing each 
product with the terms of (49) into partial fractions, one readily notices that 
soft denominators in the prefactors require like signs. Therefore, we still end 
up with exactly two contributions to the hard thermal loop, a' — a\ = a2 ~ 
±1. Moreover, of all the various combinations only those contribute to hard 
thermal loops where propagators of different statistics are combined so that 
in the end there is always a sum n + n as in (40), for the difference in the 
arguments of ft and n does not matter. 

The same argument applies to higher vertex functions, since in all the fermionic 
Green functions there is just enough power from hard loop momenta to pro- 
duce a hard thermal loop oc when all the energy denominators up to one 
are soft. 

Usually these two contributions can be identified since only the first term of 
the r.h.s. of (41) is kept and this changes sign with the spatial loop integration 
variable. For the improved hard thermal loops, however, we have to keep also 
the subsequent term in (41), which does not. 

Let us finally write down one example for an improved HTL vertex with 
external fermions. Either by direct calculation or by functional differentiation 
of (47) one finds for the quark-quark-gluon vertex 



- gT-f^Mf I (yg + ^) (yr + ^) 



+ 



{YQ - Sf) (YR - ^ 



(50) 



Just as the analytic structure of the improved fermion self-energy turned out 
to be modified in a somewhat different fashion than was the case for the gluon 
self-energy, the vertex functions involve shghtly different functions. Neverthe- 
less, these can again be expressed in terms of integrals involving the function 
A4 introduced in (25). For instance. 



^ tr (y'r'^) = gTM'f £ [m(Q+, R+) + M{Q-, R-) 



with Q+ = (go + ^,q) and g_ = (Qo 



dm 



(51) 

(52) 
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Again the inclusion of the asymptotic mass (39) smoothes out the singularities 
in A4 in the final integration over a. 



5 Conclusion 

To summarize, dressing the hard propagators of the hard thermal loops by 
the asymptotic masses oc gT that pertain to the transverse branch of the 
gluonic excitations and to the normal branch of the fermionic ones removes all 
coUinear singularities of the hard thermal loops and also preserves their gauge 
invariance. Moreover, the elegant effective actions of Braatcn and Pisarski 
and of Taylor and Wong could be generalized to summarize the improved 
hard thermal loops in an equally compact form. 

The gauge invariance of the improved hard thermal loops appears to be par- 
ticularly encouraging to use them in place of the original hard thermal loops, 
where the latter lead to singular results when external lightlike momenta are 
involved. However, we have not yet shown that the systematics of a resummed 
perturbation theory built on the now everywhere well-defined hard thermal 
loops is as it was. It may well be that the would-be coUinear singularities, 
which are only logarithmic at the level of hard thermal loops, build up in 
higher loop diagrams to a degree that overpowers the suppression by pow- 
ers of g. Indeed, in Ref. [16] it has been found that in QCD the resummed 
one-loop diagrams are sufficiently singular to contribute to the sublogarithmic 
terms in eq. (19). Actually, the imaginary part (20) is even less stable and be- 
comes large already for timelike momenta with Q"^ /q^ ~ (?, thus preventing the 
longitudinal plasmon branch from coming arbitrarily close to the light-cone. 
This behaviour, which is completely opposite to the one observed in the case 
of scalar electrodynamics [14] , will be the subject of a separate investigation 
[23]. 

Another place where coUinear singularities do not cancel from a resummed 
calculation using ordinary hard thermal loops is the case of real soft photon 
production [13]. Using the improved hard thermal loops, on the other hand, 
gives a finite result when calculating the soft contribution at resummed one- 
loop order 

^dW ^^T^m\J^]J^]. (53) 

C«3p softcontr. 277^ \EJ \g J \Mf J ' ^ ' 

which coincides with the leading logarithms of Ref. [24] . 

However, there is also a hard contribution which has to restore independence 
of the scale A separating soft from hard momenta. For this it has to be such 
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that its mass singularities are also cut off in a way that produces a In^l/g) 
besides a ln(T/A). Just like the internal hard propagators in the improved hard 
thermal loops this requires hard diagrams which are dressed by higher order 
diagrams. Again we expect that the main effect will be from the asymptotic 
thermal masses, which, as we have mentioned before, do not depend on a low- 
momentum limit (but they depend on A, which would become important in 
a more accurate calculation). So the same mechanism that renders finite the 
soft contribution should also apply to the hard one, and the last logarithm in 
(53) should combine into ln(T/M/) r^\n{l/g). 

Clearly, much work is still needed to first verify that an improved resummed 
perturbation theory really works as sketched here, and second to fully evaluate 
the sublogarithmic contributions. Encouraged above all by the gauge invari- 
ance of the improved hard thermal loops, we expect them to play a central 
role in this and analogous problems. 
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Appendix 

In order both to corroborate the results obtained in sect. 3.1 and to show a 
potential pitfall in their derivation, let us recalculate (15) by writing first 



(This assumes transversality of H, which is guaranteed by the gauge invariance 
of the (improved) HTL effective action.) 

The potential HTL contributions of U^i, give 



noo(Q) = -^n, with Ui^B^'^u 



(54) 





{PQ) 



21 



(55) 
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This can be rewritten as 





N (pq)^l 


^ V 





+ e2(4m^-g2)|?^A„A- (56) 



{I) 



{II) 



where, superficially, only the first part appears to be a hard thermal loop. It 
yields 



\{i)(n\ — ^ / (lppn(p) I dz 



W{Q) = 



g2 27r2 



2 ■ 



(57) 



The result (57) has exactly the form one obtains by expanding the manifestly 
gauge invariant effective action of eq. (26), which involves one more denomina- 
tor than the one eq. (34). However, the derivative n' implicit in the averaging 
in (26) is replaced by —2n. For conventional hard thermal loops, this does not 
make any difference. But close to the hghtcone it does. For example, the result 
(19) is not exactly reproduced — the -\-\ is missing. 

What goes wrong here is that one has to divide II^ by in order to obtain 
Hoo. This way it happens that, close to the light-cone, the second contribution 
in (56) can no longer be neglected. It reads 



2 

nf)(g) = (4m^-Q2)^ / dp 



n{p) 



p Ql{l + ml^/p^) - q^' 



(58) 



Indeed, this by itself does not give rise to a HTL contribution, except when 
the in the numerator is removed, which is done in (54). Then there is a 
linear singularity for Q"^ ^ Q that is cut off by m^^jp?' in the denominator, 
which heaves two powers of hard momentum into the numerator. This precisely 
accounts for the missing contribution to noo((5^ = 0). Moreover, adding (58) 
to (57) can be shown to be equivalent to replacing n by in the latter so 

that the gluon self-energy is exactly as prescribed by the manifestly gauge- 
invariant effective action (26). 
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